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Abstract 

Aim: Our study aimed to investigate the sulfur compounds’ antidiabetic effect in Ferula latisecta. 

Material and Methods: The molecular docking method investigated the a-glucosidase inhibitory effect of the 
components. In our study, the pharmacokinetic properties of Ferula latisecta compounds were also investi- 
gated with the SwissADME method, and the toxicity risk analyzes were investigated with Protox II tools. Ferula 
latisecta compounds were drawn from the literature in Chemdraw and and a-glucosidase enzyme structure 
was obtained from Protein Data Bank. Finally, the molecular interaction analysis between a-glucosidase and 
compounds from Ferula latisecta was performed by AutoDock 1.5.7, Molecular interactions were investigated 
using Discovery Studio Visualizer and Ligplot 2.1 program. 
Results: All the selected sulfur compounds from Ferula latisecta followed Lipinski's rules, had sufficient binding 
energy, and lacked toxicity; therefore, they were appropriate candidates for a-glucosidase inhibition. Among 
these compounds, 2-(4-hydroxyphenyl) ethyl lignocerate and isosco-poletin showed the lowest binding energy 
and the highest inhibitory effect on a-glucosidase enzyme with -91 and -77 kcal/mol, respectively. 

Discussion: These compounds also indicated a lower binding energy than the standard inhibitor (miglitol). 
Among the sulfur compounds in Ferula latisecta 2-(4-hydroxyphenyl) ethyl lignocerate and isoscopoletin were 
predicted to be the potent inhibitors due to having more hydrogen bonds and hydrophobic interactions with the 
active site of a-glucosidase. 
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Introduction 

Diabetes mellitus (DM) continues to be a significant threat to human 
health today, and it is estimated that by 2045, the number of patients 
worldwide will exceed 642 million [1]. DM is a chronic metabolic disease 
characterized by high blood sugar levels due to insulin resistance 
[1,2]. High blood glucose levels also damage blood vessels and nerves, 
causing various health problems such as hypertension, cardiovascular 
disease, blindness, stroke, amputations, kidney, and dental conditions 
[3,4.5]. Although different drug treatments are available today, the side 
effects of drugs continue to be a problem. 
Antidiabetic drugs fall into several categories, including sulfonylureas, 
bioguanidines, insulin mimetics (glucagon-like peptide analogs), and 
a-glucosidase inhibitors. 

The discovery of new drugs is essential because of the side effects 
of existing antidiabetic drugs. Academic studies have recently 
confirmed medicinal plants’ beneficial effects in treating diabetes 
[6,7]. Researched impacts of plants include improving glycemic control, 
lowering serum lipid levels, inhibiting oxidative stress, and improving 
inflammatory response [8,9]. 

The genus Ferula is mainly distributed throughout central and 
Southwest Asia (especially Iran and Afghanistan], the Far East, North 
India, and the Mediterranean. The main phytochemicals present in 
the genus are as follows: coumarin, coumarin esters, sesquiterpenes, 
sesquiterpene lactones, monoterpene, monoterpene coumarins, 
prenylated coumarins, sulfur-containing compounds, phytoestrogen, 
flavonoids, and carbohydrates [10]. 

People used these plants’ roots, leaves and fruits as vegetables, spices, 
or medicine. In addition, other species of this genus have been reported 
to be popularly used for sedatives, digestive disorders, rheumatism and 
arthritis, neurological disorders, inflammations, dysentery, headaches, 
and toothaches [11,12]. 

Ferula latisecta (F latisecta) has been used in Iranian folk medicine to 
treat parasitic diseases, relieve stomachaches in infants, and to control 
diabetes. In other studies, it has been shown that essential oil obtained 
from the above-ground parts of F. latisecta has an antimicrobial effect 
13,14]. 
na previous study, streptozotocin-induced diabetic male Wistar rats 
were treated with F. latisecta root (400 mg/kg/day) for four weeks, 
resulting in significantly reduced LDL levels in the kidney and liver (p < 
0.05) and thiol (p < 0.05) and superoxide dismutase (p < 0.01). In addition, 
the root of F. latisecta was found to reduce serum total cholesterol 
evels (p < 0.05} and prevent the progression of hyperglycemia [15]. 
The mechanism of action of many drugs used to treat diabetes 
occurs through the inhibitory effect of a-glucosidase (AG). AG is 
a carbohydrate hydrolyzing enzyme. Overactivity of a-glucosidase 
causes glucose to be absorbed by the small intestinal lumen and 
enter the bloodstream. Glucose release results in increased blood 
glucose levels and uncontrolled hyperglycemia in type-2 DM patients. 
a-Glucosidase is one of the glucosidase enzymes responsible for the 
mechanism of DM. a-glucosidase catalyzes the hydrolysis reaction by 
breaking the a-(1,4)-glycosidic bond of carbohydrates to become a free 
monosaccharide {a-D-glucose) before entering the bloodstream. The 
AG enzyme plays a vital role in carbohydrate metabolism by converting 
starch and disaccharides into glucose [16]. As a result of inhibiting the 
enzyme a-glucosidase, carbohydrate digestion slows down and causes 
glucose to enter the bloodstream later. As a result, it constitutes a 
proper prophylactic treatment for hyperglycemia [17]. Therefore, 
nhibition of AG is a practical therapeutic approach for treating DM. 
Our study aims to determine the toxicity of eight compounds [12] 
of F. latisecta by in silico methods and investigate the anti-diabetes 
effect by examining the inhibitory effect of the compounds on the 
a-glucosidase enzyme by molecular docking method. 
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Material and Methods 
First, important sulfur-containing phenolic compounds were obtained 
from the literature on F. latisecta [12]. Next, compounds were drawn 
using ChemDraw software. Finally, a suitable crystallographic 
structure with a resolution of 1.60A with the code PDB ID 3A4A from 
the a-glucosidase enzyme protein data bank was obtained from the 
protein database (https://pubchem.ncbi.nim.nih.gov/). 

Determination of pharmacological properties of compounds by 
Lipinski parameter 

The potentially effective sulfur compounds of F. latisecta were 
evaluated using the Lipinski parameter to inhibit the activity of the 
a-glucosidase enzyme. All the compounds showed the properties 
required by Lipinski's parameter; therefore, they were predicted to 
have optimal adsorption. The SwissADME database (http://www.swiss 
adme.ch/ index .php) was used to obtain the Lipinski properties of the 
compounds. 

Evaluating the toxicity of the selected compounds 

F. latisecta contains many sulfur compounds. Among the most 
highlighted sulfur compounds, latisulfide A, latisulfide B, latisulfide 
C, latisufide D, and latisulfide E, t-cadinol, 2-(4-hydroxyphenyl} ethyl 
lignocerate and isoscopoletin are also found in this plant. Phenolic 
compounds: Lack of toxicity is one of the critical factors for choosing 
a compound as a therapeutic candidate. Therefore, in this study, 
the toxicity of each of the sulfur compounds, such as liver toxicity, 
carcinogenicity, immunotoxicity, mutagenicity and the toxicity class 
of the compounds were examined using the Protox tool (https://bio. 
tools/protox} [18]. In addition, the accuracy of the toxicity analysis 
performed with the PKCSM software (https://biosig.lab.ug.edu.au/ 
pkcsm/prediction) was checked. 
Molecular interaction analysis 
This study used the Autodock vina tool (version 1.5.7) [19] to investigate 
the molecular interaction between a-glucosidase enzyme and 
selected ligands. Before the docking analysis, the enzyme structure 
was optimized by removing excess ligands and water molecules using 
the BIOVIA Discovery Studio 2021 program. The Spartan 14 (Version 
114) program optimizes all compounds for energy. Kollman charges 
were determined by adding polar hydrogens to the protein using 
the AutoDock vina 1.5.7 tool. The partial load of the compounds was 
calculated using Compute Gasteiger in AutoDock 1.5.7 tool. 
The x, y, and z coordinates were determined to bind the a-glucosidase 
enzyme to its catalytic site (x:20.53, y: -1011 2:22.38 x:40, y:-40 z-40). 
Finally, molecular interactions and binding types between the selected 
compounds and the a-glucosidase enzyme were investigated using the 
Discovery Studio visualizer and Ligplot (version 2.2.8) programs [20]. 


Results 

Pharmacological properties of the selected compounds 

The selected ligands: latisulfide A, latisulfide B, latisulfide C, latisufide 
D, latisulfide E, t-cardinal, 2-(4-hydroxyphenyl) ethyl lignoceric 
and isoscopoletin were obtained from the literature and drawn in 
ChemDraw software [12]. 

Toxicity analysis of the selected compounds 

The compounds were also evaluated for their toxicity and toxicity 
class. The criterion for selecting a compound as a drug candidate is 
the safety of the compound, and the compound should not show any 
toxicity. 

Compounds contained in F. latisecta were evaluated in terms of 
toxicity and toxicity class. In the ADME study conducted on SwissADME 
software (http://www.SwissADME.ch/], it was found to be appropriate 
among the compounds in terms of pharmacokinetic properties, except 
for 2-(4-hydroxyphenyl} ethyl lignocerate. Toxicity estimation was 
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evaluated with the protox II tool program. According to the protox II tool 

prediction program, it was observed that hepatotoxic and immunotoxic ‘ : ; 

effects may occur in all compounds. However, when we look at the j a 
probability values, it has been determined that the risk of toxicity is = 
low. Toxicity class 4 indicates low toxicity of the ingredients. When 4 
the PKCSM program, which is a different software, was examined ——— A | 
separately, no toxicity was detected. 


Molecular interaction analysis 
The results of the molecular interaction between the selected compound 


Table 1. Investigation of pharmacologic parameters (Drug-like) for = 
sulfur compounds in F. latisecta according to Lipinski rule. & 
ais a = a 
Rea Hydrogen bond Hydrogen bond Molecular (<5) — - 
P donors (<5) acceptors (s10) mass (<500) g afn 
a m iii | 

Latisulfide A 0 4 220.36 3.24 = = 

Latisulfide B 0 4 246.40 3.80 - = S 

Latisulfide C 0 4 26042 419 —_—a ——— 

— — ‘ hy 
Latisulfide D 1 3 178.32 2.67 
d 

Latisulfide E 1 4 35449 4.65 7 

t-cadino! ‘I 1 222.37 3.18 

2 (etycronphey 1 3 460.74 9.30 Figure 1. Two-dimensional (2D) and three-dimensional (3D) plots of 

on F i = ape a-glucosidase interactions with selected potential inhibitors. 


a) 2D representation of the molecular interactions between enzyme 
Table 2. Evaluation of toxicity and toxicity class of the selected Ferula  2-(4-hydroxyphenyl) ethyl lignocerate b) 3D representation of 
latisecta compounds. molecular interactions between enzyme and 2-(4-hydroxyphenyl) 
ethyl lignocerate. c) 2D representation of molecular interactions 
between enzyme and isoscopoletin. d) 3D representation of molecular 
interactions between enzyme and isoscopoletin. 


Compound Hepatotoxicity Carcinogenicity Immunotoxicity Mutagenicity 


Latisulfide A H NC NM 
Latisulfide B H NC NM 
Latisulfide C H NC NM 
Latisulfide D H NC NM 
Latisulfide E H NC NM : 
t-cadinol H NC NM 
2-(4-hydroxyphenyl) 7 
ethyl lignocerate 7 NM il 
H: hepatotoxicity, NC: no carcinogenicity, NM: no mutagen, |: immunotoxicity 
Table 3. Interaction and binding energy of sulfur compounds with 
amino acids of a-glucosidase. 
Binding energy Hydrogen Hydrophobic 
Compounds , , : 5 
P (Kcal/mole) interaction interaction 
Latisulfide A -6.2 Arg442 Phe303, His280 
Aes Phe178, Arg315, 
Latisulfide B -6.5 - His280 
pars 7 Lys156, Arg315, Tyr158, 
Latisulfide C 6.6 Arg315 Phet59, Phet78 
aoe: yri58, Phe1/8, 
Latisulfide D -5.5 Asp215, Asp352 val216, Phe159 
Latisulfide E 7 Tyr407, Glu405 Trp36 a 
t-cadino “11 - Spal Wiel ai 


Phe3t4, Arg315 


2-(4-hydroxypheny!) 


ethyl ignocerate 91 Pro66 Protst, Trp238, Vald10 Figure 2. Three-dimensional representation of the receptor and 
5er236, Asn236, 6lu429 isoscopoletin ligand surface interactions in the active site of 
isoscopoletin “iT “Asn3i7 t le419, His423 5 - : 2 4 
a-glucosidase regarding polarity, hydrogen bonding, hydrophobic and 
Miglitol (Standard compound) Al GIn353, Arg442 ionized surface of the exchanges 
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and the a-glucosidase enzyme obtained by molecular coupling assay 
are given in Table 1. All compounds showed good binding with the 
enzyme. The binding energy, hydrogen and hydrophobic interactions 
between the compounds and the enzyme are also presented in Table 3. 
The results of the molecular docking study showed that all sulfur 
compounds in F. latisecta can bind to the active site of a-glucosidase 
and inhibit the activity of this enzyme. According to the docking analysis, 
the binding energies of the studied compounds are different; therefore, 
the critical energy ranges from -5.5 to -91 kcal/mol. The lower the 
binding energy level (negative), the stronger the binding between the 
receptor (enzyme) and ligands (compound or inhibitor). Among the 
selected sulfur compounds, 2-(4-hydroxyphenyl) ethyl lignoceric with 
a value of -91 kcal/mol and isoscopoletin with a weight of -77 kcal/mol 
showed the lowest binding energy. Therefore, they are estimated to 
provide the highest inhibitory effect among the compounds. 

Previous studies have shown that amino acids such as LYS-156, GLN- 
279, HIS-280, PRO-312 GLU-277, ASP-352, ASN-415, and ARG-442 play a 
key role in the interaction of enzyme and inhibitor in the active site 
of a-glucosidase [21]. Miglitol was used as a standard a-glucosidase 
inhibitor. The insertion results showed that miglitol, with a binding 
energy of -71 kcal/mol, has a significant binding energy and hydrogen 
bonding with amino acids GIn353 and Arg442. 

Pro66 (3.60 and 5.50 A) 2-(4-hydroxyphenyl) ethyl lignoceric complex 
and isoscopoletin complex with amino acids Ser236, Asn235, Glu429, 
Asn317 also have hydrogen bonds (Figure 1 a,b). Isoscopoletin showed 
more hydrogen bond interactions than other compounds. These two 
compounds have lower binding energies than other compounds. This 
also shows that it can bind more strongly to the enzyme. Figure 1 shows 
the three-dimensional view of the ionization, hydrophobic, polar, and 
hydrogen bonding reactions between a-glucosidase and isoscopoletin. 


Figure 3. Comparison of isoscopoletin and 2-(4-hydroxyphenyl) 
ethyl lignocerate miglitol (standard inhibitor) in the active site of 
a-glucosidase enzyme. 
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Discussion 

In the present study, 2-(4-hydroxyphenyl) ethyl lignoceric has shown 
that F. latisecta has the lowest binding energy among the phenolic 
compounds. In another study, while the synthesis and _ inhibitory 
potentials of sulfur compounds were investigated, it was found that 
compounds such as methazolamide, acetazolamide, and timolol 
(having sulfur atoms with binding energies of -4.3, -4.8 and 4.5, kcal/ 
mol in their structure, respectively) inhibited the a-glucosidase 
enzyme significantly [22]. In this study, with a -91 kcal/mol value, 
2-(4-hydroxyphenyl} ethyl lignocerate showed the best inhibition value 
among sulfur compounds. 

In addition, miglitol did not show any hydrophobic interaction with 
a-glucosidase. However, it showed hydrophobic interactions with 
2-(4-hydroxyphenyl} ethyl lignoceric and a-glucosidase Prot51, Trp238, 
and Val410. Exchange of isoscopoletin with a-glucosidase showed 
hydrophobic interaction with amino acids lle419 and His423. 
Conclusion 

When we looked at the results of our study, it was clear that the 
pharmacokinetic properties of important sulfur compounds were 
good. Still, hepatotoxicity and immunotoxicity risk were observed in all 
components; however, the low toxicity probability levels and the toxicity 
class 4 show that it does not carry a significant stake in toxicity. Also, 
the hepatotoxic effects of hepatotoxic compounds were examined in a 
different software PKCSM program. When the data obtained from this 
program were analyzed, it was seen that they were not toxic. Based 
on the results of the docking study, it can be concluded that of the 
eight compounds in F. latisecta, 2-(4-hydroxyphenyl) ethyl lignocerate 
and isoscopoletin are more effective in inhibiting the a-glucosidase 
enzyme. These compounds also exhibited lower binding energies 
compared to miglitol. Therefore, they are likely to be more potent 
inhibitors than miglitol. Considering that in silico studies are used for 
preliminary studies and predictions, it is thought that the results of our 
study should be validated with in vivo and in vitro studies. 
Acomparison of 2-(4-hydroxyphenyl) ethyl! lignocerate and isoscopoletin 
with miglitol is shown. Both selected compounds appear to interact 
with different amino acids in the catalytic pocket of the a-glucosidase 
enzyme. 


Scientific Responsibility Statement 

The authors declare that they are responsible for the article's scientific content including 
study design, data collection, analysis and interpretation, writing, some of the main line, or 
all of the preparation and scientific review of the contents and approval of the final version 
of the article. 


Animal and Human Rights Statement 

All procedures performed in this study were in accordance with the ethical standards of the 
institutional and/or national research committee and with the 1964 Helsinki Declaration and 
its later amendments or comparable ethical standards. 


Funding: None 


Conflict of Interest 
None of the authors received any type of financial support that could be considered 
potential conflict of interest regarding the manuscript or its submission. 


References 

1. Cho NH, Shaw JE, Karuranga S, Huang Y, da Rocha Fernandes JD, Ohlrogge AW, et al. IDF 
Diabetes Atlas: Global estimates of diabetes prevalence for 2017 and projections for 2045. 
Diabetes Res Clin Pract. 2018:138:271-81. 

2. Association AD. 7. Diabetes Technology: Standards of Medical Care in Diabetes—2021. 
Diabetes Care. 2020:44(Suppl. 1}:S85-99, 

3. Cui Y, Zhang L, Zhang M, Yang X, Zhang L, Kuang J, et al. Prevalence and causes of low 
vision and blindness in a Chinese population with type 2 diabetes: the Dongguan Eye Study. 
Sci Rep. 2017:7(1):11195. 
4. Helgason CM. Blood Glucose and Stroke. Current Treatment Options in Cardiovascular 
Medicine. 2012;14(3):284-7. 
5. Zhang P. Gregg E. Global economic burden of diabetes and its implications. The Lancet 
Diabetes and Endocrinology. 2017;5(6):404-5. 


Molecular interactions of Ferula latisecta 


6. Ghorbani A. Best herbs for managing diabetes: a review of clinical 
studies. Brazilian Journal of Pharmaceutical Sciences. — 2013;49(3):413-22. 
7. Hosseini A, Shafiee-Nick R, Ghorbani A. Pancreatic beta cell protection/regeneration with 
phytotherapy. Brazilian Journal of Pharmaceutical Sciences. 2015;51(1): 1-16. 

8. Ghorbani A. Mechanisms of antidiabetic effects of flavonoid rutin. Biomed Pharmacother. 
2017; 96:305-12. 

9. Prabhakar PK, Doble M. Mechanism of action of natural products used in the treatment of 
diabetes mellitus. Chin J Integr Med. 2011;17(8):563-74. 

0. Sonigra P Meena M. Metabolic Profile, Bioactivities, and Variations in the Chemical 
Constituents of Essential Oils of the Ferula Genus [Apiaceae]. Front Pharmacol. 
2021:11:608649. 

‘1. Dehghan G, Shafiee A, Ghahremani MH, Ardestani SK, Abdollahi M. Antioxidant Potential 
of Various Extracts from Ferula szovitsiana. in Relation to Their Phenolic Content. 
Pharmaceutical Biology. 2007;45(9)-691-9. 

2. Soltani S, Amin GR, Salehi-Sourmaghi MH, Schneider B, Lorenz S, lranshahi M. Sulfur- 
containing compounds from the roots of Ferula latisecta and their cytotoxic activities. 
Fitoterapia. 2018;124:108-12. 

3. Habibi Z, Salehi P. Yousefi M, Hejazi Y, Laleh A, Mozaffarian \V, et al. Chemical composition 
and antimicrobial activity of the essential oils of Ferula latisecta and Mozaffariania insignis 
from Iran. Chemistry of Natural Compounds.2006;42(6):689-92. 

4. Iranshahi M, Hassanzadeh-Khayat M, Bazzaz BSF, Sabeti Z, Enayati F High Content of 
Polysulphides in the Volatile Oil of Ferula latisecta Rech. F et Aell. Fruits and Antimicrobial 
Activity of the Oil. Journal of Essential Oil Research.2008;20(2):183-5. 

5. Javanshir S, Soukhtanloo M, Jalili-Nik M, Yazdi AJ, Amiri MS, Ghorbani A. Evaluation 
Potential Antidiabetic Effects of Ferula latisecta in Streptozotocin-Induced Diabetic Rats. J 
Pharmacopuncture. 2020;23(3):158-64. 

6. Abdjan MI, Aminah NS, Kristanti AN, Siswanto |, Ilham B, Wardana AP. et al. Structure- 
based approach: molecular insight of pyranocumarins against a-glucosidase through 
computational studies. RSC Advances. 2023:13(6):3438-47. 

7. Gong L, Feng D, Wang T, Ren Y, Liu Y, Wang J. Inhibitors of a-amylase and a-glucosidase: 
Potential linkage for whole cereal foods on prevention of hyperglycemia. Food Sci Nutr. 
2020;8(12):6320-37. 
8. Drwal MN, Banerjee P Dunkel M, Wettig MR, Preissner R. ProTox: a web server for the in 
silico prediction of rodent oral toxicity. Nucleic Acids Res. 2014;42(Web Server issue]:W53-8. 
9. Trott 0, Olson AJ. AutoDock Vina: improving the speed and accuracy of docking with a new 
scoring function, efficient optimization, and multithreading. }) Comput Chem. 2010;31(2):455- 
61. 
20. Laskowski RA, Swindells MB. LigPlot+: Multiple Ligand-Protein Interaction Diagrams for 
Drug Discovery. | ChemInf Model. 2011;51(10):2778-86. 

21. Zhang J, Li YN, Guo LB, He J, Liu PH, Tian HY, et al. Diverse gallotannins with a-glucosidase 
and a-amylase inhibitory activity from the roots of Euphorbia fischeriana steud. 
Phytochemistry. 2022;202:113304. 

22. Gollapalli M, Taha M, Javid MT, Almandil NB, Rahim F, Wadood A, et al. Synthesis of 
benzothiazole derivatives as a potent a-glucosidase inhibitor. Bioorg Chem. 2019;85:33-48. 


Eurasian Clinical and Analytical Medicine 


2 


